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Abstract 



In this work we analyzed the interaction of isolated solitary structures and 

ion-acoustic radiation. If the radiation amplitude is small solitary structures 

persists, but when the amplitude grows energy transfer towards small spatial 

scales occurs. We show that transfer is particularly fast when a fixed point of 

a low dimensional model is destroyed. 
05.45.+b 
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I. INTRODUCTION 



Langmuir turbulence has been one of the most studied problems in modern nonlinear plasma 
physics. Over the last years a great deal of effort has been directed to its analysis, as well as 
to the analysis of related subjects as soliton dynamics, collapse, nucleation of cavitons, elec- 
tromagnetic emission, and others PJ. More recently, attempts have been made to understand 
the turbulence in terms of concepts of nonlinear dynamics and chaos &-§|. 

The conservative version of Langmuir turbulence is described by the Zakharov equations 
that couple the slowly varying amplitude of a high-frequency electric field, the Langmuir 
field, to slow density fluctuations, the ion-acoustic field. Decay processes deposit energy into 
Langmuir fluctuations with long wavelengths and if the energy thus accumulated exceeds 
the threshold for modulational instability, solitons can be formed. 

In addition to solitons a certain amount of ion-acoustic radiation is also generated, a 
fact that creates the possibility of nonlinear wave interaction involving these two types 
of structures: solitons of the Langmuir field and long wavelength ion-acoustic radiation. 
In more specific terms, what happens is that as solitons are formed their shapes exhibit 
temporal oscillations 0; if ion-acoustic fluctuations are also present, the possibility exists 
of interaction between the oscillatory degrees-of-freedom of solitons and the oscillating ion- 
acoustic waves. One has two length scales in the region of long scales. One of them is the 
soliton length scale, we shall call it L s , and the other is the length scale of the ion-acoustic 
fluctuations, Lj. It has been shown that depending on the general conditions of the system, 
the mentioned interaction may lead to intense energy transfer from the spectral region of long 
wavelengths to the region with much shorter wavelengths, we call it L s h with L s h <C L s , Lj. 
As energy moves into modes with small wavelengths, dissipation becomes progressively more 
important. However, as we are interested only in nonlinear transfer processes, we discard 
dissipation in a first approximation. It has been suggested that energy transfer occurs when 
the interaction is of chaotic nature. Presumably the process underlying the transfer is related 
to the diffusive processes induced by the presence of a stochastic drive in the system; the 
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stochastic drive would be formed by the chaotic degrees-of- freedom ||. 

Now we come to our point. In various earlier simulations ||,[7|,|§ a modulationally per- 
turbed plane wave is launched into the system. If the system is unstable a number of solitons 
and additional ion-acoustic radiation are formed. Solitons interact with each other and with 
the radiation, and transfer of energy towards small spatial scales L s h may take place if 
nonintegrable features are prominent. The problem here is that this kind of simulations 
does not examine properly the interaction of individual solitons and the radiation, since 
soliton-soliton collisional processes cannot be disregarded under such conditions. It is not 
even clear which type of interaction, if soliton-soliton or soliton-radiation, is the dominant 
responsible for the transfer. In addition, several systems display a small soliton density that 
collisions are unlikely - in these systems one should focus attention on the individual inter- 
action involving one single oscillating soliton and ion-acoustic waves. This is the purpose 
of the present paper. We shall examine the system evolving from an initial condition where 
only one single oscillating soliton and some radiation are present. In addition to simulations 
we develop a model where we perform averages over fast variables in order to make estimates 
with regard to the behavior of the collective variables of the system. 

As will become clear, energy transfer starts to take place when the collective variables 
become chaotic. In general we shall see that while for moderately and small amplitudes of 
the perturbing ion-acoustic radiation solitons can be at least seen as metastable structures 
in the system, for large amplitudes transfer is fast and initial solitary structures are rapidly 
destroyed. It has been argued that soliton are robust enough to describe final states of this 



type of system ||lTf-|T3|j . But what we see here is that even if some solitons are present in 
asymptotic states of large amplitude regimes, those solitons are not the same present in 
earlier times - after the initial solitons are destroyed there are long stretches of times over 
which no organized structures are seen. 

We finally mention that a number of works have already analyzed the interaction of 
localized structures and perturbations with longer wavelengths. In some of them only low 



dimensional models were investigated |L4]], and in others where full simulations were per- 



formed, chaotic dynamics was not the issue, although some nonintegrable features like soliton 
fusion have been reported ||. 

We organize the paper as follows: in §2 we introduce the basic model and the numerical 
techniques to be used here; in §3 we discuss our initial conditions and perform the appropriate 
averages to single out the relevant collective variables; in §4 we compare the low dimensional 
model with full simulations, and in §5 we summarize the work. 

II. BASIC EQUATIONS AND NUMERICAL TECHNIQUES 

The one dimensional Zakharov equations governing the Langmuir turbulence can be written 
in the adimensional form || 

id t E + d 2 x E = n E, (1) 

d 2 t n-dln = dl\E\\ (2) 

with dt = d/dt, d x = d/dx. E(x,t) is the slowly varying amplitude of the high-frequency 
Langmuir field and n(x,t) are slow density fluctuations associated with the ion-acoustic 
field. The NLS equation 

id t E + d\E + |£| 2 E = (3) 

is obtained from the set (|1|), (0) if one is allowed to approximate Eq. (0) in order to replace 
n with — l-E 1 ] 2 + const. This approximation is called subsonic because it requires very slow 
time scales that dtn(x,t) -C 8%.n(x,t). 

Our numerical approach is based on a pseudo-spectral method. We assume spatial peri- 
odicity with basic length L and expand E(x,t) and n(x,t) into Fourier series as 

+— 

E{x,t)= E E m {t)e imkx , and 
m =-f 

+f 

n(x,t)= E n m {t)e imkx . (4) 



The basic wavevector is defined in terms of the system length L as k = 2n/L, and the integer 
N represents the number of modes used in the simulations. To represent a continuous system 
one should take the limit iV — > oo. In practice we let iV = 1024, removing half of the modes 
to cure aliasing problems associated with the FFT routines. Comparisons with iV = 2048 
indicates numerical convergence in terms of number of modes. Accuracy is further checked 
by varying the tolerance factor of the numerical integrator and by monitoring the conserved 
energy [[j]. We find that relative fluctuations in energy are about one part in 10 6 — 10 8 and 
that variations of tolerance factor do not produce alterations in the outcome of runs. 

Solitons of amplitude v2 a s are formed when an homogeneous train of Langmuir radiation 
of amplitude \E Q \ 2 ~ a s becomes modulationally unstable. The subsonic growth rate V for a 
perturbation with wavevector k superimposed on the homogeneous train can be estimated 
as 

T - ~ ^\E Q \ 2 -k 2 . (5) 

From relation (|5|) one sees that the only unstable modes are those for which \E \ 2 > k 2 . Now 
when \E \ 2 — k 2 <C 1, T <C k. If this condition holds for the majority of modes, d t <C d x , 
ion-acoustic fluctuations are mostly enslaved to the Langmuir field, and approximation ([|) 
can be used. On the other hand when \E \ 2 is not exceedingly small ion-acoustic fluctuations 
with k -C \E \ may not be completely enslaved to the Langmuir field. Those free fluctuations 
are to be seen as independent degrees-of-freedom whose presence is capable of destroying 
the integrability of the system. Given that the maximum growth rate occurs for k max ~ \E Q \ 
and that the typical length scale of a soliton is given by L s ~ l/k max , free ion-acoustic 
radiation of wavevector fc* = 1/Lj typically appears in the spectral region for which 

U > L s . (6) 
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III. COLLECTIVE VARIABLES AND LOW DIMENSIONAL MODEL 



Our system is multidimensional but we would like to see whether a small subgroup of modes 
is more active than the remaining. If this is the case one could try to describe the basic 
features of the full dynamics by a low dimensional approximation. As it turns out, such an 
approximation appears to be possible. 

To see how to obtain the low dimensional model, we proceed as follows. We first recall 
that as initial conditions we are interested in configurations with isolated solitary structures. 
To represent this sort of states either analytically or in the simulations we shall first deter- 
mine the stationary one-soliton-solution for the full problem. We start by taking d t = in 
Eqs. (0) and (0) from which we first get n m — \E S \ 2 + constant. Substituting this relation 
into Eq. ([!]), after some algebra one obtains 



which is the expression we are looking for. £ is an arbitrary factor that measures either the 
amplitude or the inverse width of the soliton. We point out that due to the nonlinearities 
and dispersion of the problem a precise balance between amplitude and width is needed. If 
we call a s = £ and w s = |, it is indeed seen that the following relation holds: 

a s = — . (8) 

w s 

We had mentioned that our interest is to see what could happen with the soliton when 
it starts to interact with free ion-acoustic radiation. Based on several results one knows 
already that the basic soliton solution must be allowed to display temporal oscillations. 
The problem now is how to describe those oscillations in a compact way. And the answer 
is known: one first write down an ansatz solution for the soliton field where amplitude is 
however not correlated to the width according to the static relation (J|). The ansatz solution 
is therefore generically written in the form 




(7) 




x 



(9) 
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where a(t), w(t), and $(t) are all unknown as yet. The phase factor $ is included to 
incorporate the complex structure of the solutions of the set ([]]), (@). As for the ion- acoustic 
field interacting with the soliton field, one writes 



n(x,t) = -\E '(x, t)\ 2 + (A(t) e lkx + c.c.) (10) 

Here we write the ion-acoustic field as a sum of the pure adiabatic response to the soliton 
field, plus some free radiation that will actually interact with the isolated nonlinear struc- 
ture. A(t) is the amplitude of the radiation field and c.c. stands for complex conjugate. 
The next step is to derive the appropriate governing equations for the four time dependent 
parameters, a(t),w(t),§(t), and A(t). This is more easily done with help of average La- 
grangean techniques. The full Lagrangean from which one obtains the original set (([]), (@) 
reads 

L = Jcdx = J [^(E*d t E - Ed t E*) - \d x E\ 2 - \E\ 2 d x v + ±[(d t v) 2 - (d x u) 2 ]]dx, (11) 

where the dynamical variable u(x,t) is introduced in the form n(x,t) = d x v(x,t). The 
Euler-Lagrange equation for E(x,t), for instance, is written as 

dC DC dC 

l d{d t E) BE x d{d x EY 1 } 

with similar expressions holding for the other variables. From expression (0) one obtains 
the complex conjugate of Eq. (p. In terms of averaged Lagrangeans, what has to be done 
now is to substitute into Eq. ([□]) the one-soliton solution, Eq. (^|), plus the ion-acoustic 
field, Eq. (|T0"D. Doing this and performing the spatial integrations one arrives at 



AW 2 2W W 2 w 2 ■ 
+ 0.429 3.290W ww A 



nA 2 tiA 2 . . 

+ ^~— (13) 



3w 3w 2 ' w 

with 1] = a(t) 2 w(t). The various numerical factors appear in Eq. (|T3| ) as a result of the 
integrals involving trigonometric and hyperbolic functions. 

Euler-Lagrange equation with respect to the variable $ indicates that r] is a constant 
of motion. As a matter of fact this feature has been already used to simplify the form of 



the Lagrangean (|T3| ) by dropping terms proportional to f] up to positive powers. Euler- 
Lagrange variational equations are then applied to the independent variables w(t) and A(t) 
to produce a two-degrees-of-freedom conservative dynamical system. If we set A — > we 
have solutions corresponding to free oscillations of the soliton shape. One can construct a 
convenient phase-space to visualize those oscillations. This is done in Fig. (|I]) where we plot 
w(t) versus w(t). The central fixed point of the figure is simply the static soliton solution 
analytically represented by Eq. (|7|), and the curves surrounding the fixed point represent 
oscillatory modes of the soliton, each mode labeled by a particular constant energy that can 
be canonically evaluated from Lagrangean ( |i~3D with A = 0. In the absence of ion-acoustic 
free fluctuations, one can estimate the position of the fixed point, 

a s = l/w s = r], (14) 

and the oscillatory frequency around the fixed point, 



= V llg ~ fl - (15) 

Given that L s = w s = l/a s , one has L s ~ 1/uj s , and given that 2ir/uJi = Li ^> L s one 
obtains a relationship involving the frequencies of soliton and ion-acoustic waves: 

LOi < LO s . (16) 

In other words, the components of the ion-acoustic field most weakly enslaved to the Lang- 
muir field are those for which both length and time scales are much longer than the scales 
corresponding to the solitons. One shall also mention that in addition to trapped orbits 
around the fixed point, open orbits are also possible. Those would represent decaying soli- 
tons for which w —>■ oo asymptotically. The fact that one has trapped and untrapped orbits 
implies that a separatrix does exist in which vicinity some amount of chaotic activity may 
be displayed if the system is in fact nonintegrable. The role of chaos, if chaos is indeed 
present, shall be better explored in the next section. 
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IV. FULL SIMULATIONS VERSUS THE LOW DIMENSIONAL MODEL 



At this point we make use of the numerical techniques discussed in §2 to compare results 
of full one dimensional simulations with the low dimensional model developed in the pre- 
ceding section. Our full simulations give an account of the behavior of a stationary soliton 
submitted to the action of long wavelength ion-acoustic perturbations. Our purpose is to 
test the robustness of the soliton solution and see what happens when it looses stability 
due to the ion-acoustic radiation. Before embarking into the simulations it is perhaps con- 
venient to preview the basic system behavior, based on possible results obtained with the 
low dimensional model. If the parameters of the low dimensional model are such that the 
corresponding nested orbits on the phase-plane w, w are mostly regular, one can expect a 
negligible influence exerted by the ion-acoustic field on the solitary structure in the full sim- 
ulations, be this structure oscillatory or not. On the other hand it may well happen that the 
low dimensional system be nonintegrable. Should this be the case, and if low dimensional 
chaos is indeed well developed, the influence of ion-acoustic waves may be strong enough to 
destroy the solitary structure. In this case our low dimensional description may be expected 
to cease furnishing reasonable results. What is likely to happen then is that the chaotic 
low dimensional degrees-of-freedom start to act like a random drive, continuously delivering 
energy in a diffusive way to all the other dynamical modes ||. Then one may anticipate 
the soliton to decrease in intensity as its energy flows away. In addition, short wavelength 
modes are expected to grow and appear in the spectrum. We shall investigate some details 
of the transfer next. 



A. Low dimensional analysis 

Let us first explore the integrability properties of the low dimensional approximation. To 
do that we examine the surface of section obtained when we record the pair of variables w, w 
each time A = with A > 0. In the context of the low dimensional analysis we examine the 
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system as an ion-acoustic wave of initial amplitude A a is added to the central fixed point 
of Fig. (|1|) - the remaining initial amplitudes corresponding to other orbits are obtained 
with help of the condition of constant energy. This constant energy is to be obtained from 
Lagrangean ([13]). Parameters are specified in the legend of Fig. @. In both low dimensional 
and full simulation we work with a soliton of a s = ^/o.l/2 and with a perturbing ion-acoustic 
wavevector ki = 0.0257. For those parameters, 

U ~ 40L S , (17) 

and 

u s = 10.8wi. (18) 

One thus has uj s 3> u>i and L s <C Li as required by the assumptions on time and length 
scales. Relation flT8| ) in particular says that if the system is indeed nonintegrable, a period 
one island is likely to appear close to the central fixed point in the w,w phase-space. 

Examining the phase plots of Fig. (0) one sees that for small amplitudes the phase-space 
is mostly regular. However for larger values two features become noticeable: (i) the dynamics 
is indeed nonintegrable, and (ii) for large enough values of the amplitude, chaotic dynamics 
is dominant. In addition, for sufficiently large amplitudes (in the present case A Q ~ 0.145) 
the central fixed point undergoes an inverse tangent bifurcation and disappear along with 
the unstable fixed point of the period one island seen in Fig. (0b). All these features strongly 
suggest that the stochastic drive mechanism may be operative causing energy transfer into 
small spatial scales for moderately large values of the perturbation. This type of behavior 
is found for other choices of the ratios u>i/u s and Li/L s as long as relations (§) and (0) are 
respected. 

B. Full simulations 

The results of full simulations can be found in Fig. where we make plots of the space- 
time history of the field \E(x,t)\ 2 , and of the average number of active modes versus time. 
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The average number of modes is an auxiliary tool that can help to study details of energy 
transfer that are not particularly apparent in the space-time plots. The average number 



of modes is denoted by J < iV| i > for Langmuir and ion-acoustic fields, respectively, and 



defined according to the following [15 
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The primes in definition (p20|) mean that the ion modes into which energy is initially placed 
are to be excluded from the summation. We do this simply to obtain clearer results. The 
problem is that since all the initial ion-acoustic energy goes into one single mode, the statis- 
tics becomes poor if we do not make the exclusion. No problems of that sort occur with the 
Langmuir field, as solitons already involve a statistically good number of modes. 

We launch a solitary structure of shape given by Eq. (|9]), with a slight mismatch between 
a(t = 0) and l/w(t = 0) such that the soliton can oscillate initially: we choose a(t = 



0) = JO. 1/2 and w(t = 0) = 1.2/a(t = 0). For small enough values of the ion-acoustic 
perturbation Fig. (|3|a) shows that the solitary structure maintains its original amplitude 
without noticeable damping. It is seen from Fig. (§b) that for this perturbing amplitude 
the number of modes involved in the dynamics does not change significantly as time evolves. 
It should be noticed that the present modes are those used to construct the solitary structure. 

For larger values of the ion-acoustic amplitude, as in Fig. (|3]c), the soliton gradually 
damps away as time advances. Fig. (|3|d) shows that energy diffusion is now present, and that 
in the ion-acoustic field it is considerably much faster than in the Langmuir field. Diffusion 
in the Langmuir field becomes in fact almost imperceptible for even smaller perturbing 
amplitudes as we shall see a little later. Note that plateaus in the plots are formed when 
all the modes used in the simulation become involved in the dynamics - at this stage energy 
would be dissipated if we had added dissipation terms for large values of wavevectors. In 
the present case represented in Figs. (||c) and (||d), the central fixed point is still present in 
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low dimensional phase-plots, as indicated by Fig. (@b). One can therefore think in terms of 
stochastic drive models to describe this type of regime ||. Although the orbits are chaotic, 
the presence of the central fixed point offers some resistance against rapid destruction of 
the low dimensional chaotic system. This low dimensional system might therefore last long 
enough to serve as a drive delivering energy to short wavelengths modes. 

Now, if A is large enough, the localized solitary structure is rapidly destroyed as indicated 
in Figs. (||e) and Figs. (||]f). Energy is transferred to short wavelengths over short periods 
of time. We point out that this fast process occurs for ion fields intense enough to destroy 
the central fixed point of the low dimensional system, as indicated in Fig. (@c). In addition, 
diffusive time scales for both Langmuir and ion-acoustic fields become similar in this fast 
regime. Under such conditions the stochastic drive may not be a very appropriate concept 
since the life time of the solitary structure is too short. 

We emphasize, therefore, that three distinct regimes appear to be present: 

(i) If the initial perturbation is small, typically A <C 0.1, there is no diffusion whatsoever 
towards small length scales. 

(ii) For larger values of the perturbation, A ~ 0.1, diffusion is observed in both Langmuir 
and ion fields. But diffusion in the ion field is much faster. If one diminishes not too much 
the perturbing amplitude and reduces the observation time, diffusion in the Langmuir field 
becomes almost imperceptible although diffusion in the ion field can still be observed. This 
is what can be seen in Fig. (^;) where one considers a perturbing amplitude smaller, but 
of the same order of magnitude, than the one used in Fig. (||d). In general, within this 
range of perturbing amplitudes, the central fixed point of the low dimensional phase-space 
is still present. This could explain the persistence of the solitary structure seen in Fig. (0c). 
Since solitons are persistent and typically chaotic here, this regime is perhaps the most 
appropriately described by the stochastic drive. The oscillating low dimensional subsystem 
formed by the soliton and the ion-acoustic wave would excite the remaining modes of the 
system. As mentioned, diffusion is very asymmetric, being much faster in the ion-acoustic 
field. But on examining Eq. @, it is not unreasonable to say that the Langmuir field 
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term, appearing in the form d%\E(x, t)\ 2 on the right-hand side, can act similarly to a source 
delivering energy to the ion field on the left-hand side. The source-like behavior would 
enhance diffusion in the ion-acoustic field. 

(iii) Finally, when the amplitude attains sufficiently large values, A > 0.1, fast diffusion 
takes place in both fields. In contrast to the preceding case, here the time scales for diffusion 
in both fields are similar. We point out that the central fixed point of the phase-plots no 
longer exists for this range of relatively large perturbing amplitudes. Again, this could 
explain the short life of the solitary structures, as seen in Fig. (He). 



V. CONCLUDING REMARKS 

In this paper we examined the interaction of an ion-acoustic harmonic mode with a 
solitary wave of the Zakharov equations. Here the interest is to see how far can a solitary 
wave resist before it is destroyed by long wavelength radiation and how this destruction 
takes place. Although some recent works show that solitons can be stable structures even 
in nonintegrable environments |10|--|T3[| , what we see here is that if chaos is strong enough 



in the low dimensional approximations, solitons are in fact destroyed and energy transfer 
towards small spatial scales takes place. 

We have observed that the dynamics can be divided into three categories as a function 
of the amplitude of the initial ion-acoustic wave. Considering \/2a s = y/0.1, for small am- 
plitudes, A <C 0.1, there is energy transfer neither in Langmuir nor ion-acoustic fields. For 
moderately large amplitudes, A ~ 0.1, diffusion is observed mostly in the ion-acoustic field, 
and for sufficiently large amplitudes, A > 0.1, diffusion is fast and equally present in both 
fields. In the intermediary regime one can think in terms of a stochastic drive delivering 
energy to modes with short wavelengths. The drive would be formed as a result of the 
chaotic, but persistent, low dimensional dynamics. Persistence follows because for not too 
large amplitudes the central fixed point of the low dimensional phase plots is still unaf- 
fected by the interaction, which means that solitons last long enough to serve as stochastic 
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drives. While soliton turbulence may well describe the regime of intermediary amplitudes, 
it may not be quite appropriate to describe the regime of large perturbing amplitudes since 
solitons readily damp away there. The characteristics of the stochastic drive are not easy 
to be obtained because the dynamics on the chaotic space w(t),w(t) is not pendulum-like. 
Therefore, some known results on pendulum-like settings |J cannot be directly used here. 
Details and comparisons with the full simulations are currently under study. 

Recalling our initial question in this paper, the conclusion is that the interaction of 
isolated solitons and ion-acoustic radiation alone is capable of driving energy transfer. 
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FIGURES 



FIG. 1. Contour levels for the unperturbed dynamics A — > 0; 77 = y/0.1/2. 

FIG. 2. Poincare plots (w, w) of the low dimensional model with ki = 0.0257 and 77 = y/0.1/2. 
A = in (a), 0.14 in (b), and 0.16 in (c). 

FIG. 3. I (= |£(x,r;)| 2 ) and V< A^ 2 > from full simulations with ki = 0.0257 and 
a(t = 0) = y/0J/2, w(t = 0) = 1.2/o(* = 0). A = 0.005 in (a) and (b), 0.1 in (c) and (d), 
0.2 in (e) and (f), and 0.05 in (g). Time has been normalized by a factor of 5000/100 = 50 and 
space by a factor of 1024/32 = 32. 
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Fig.1 :Erichsen,Oliveira, and Rizzato 
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Fig2a: Oliveira, Erichsen, and Rizzato 
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Fig2b: Oliveira, Erichsen, and Rizzato 
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Fig2c: Oliveira, Erichsen, and Rizzato 
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Fig. 3b: Erichsen, Oliveira, and Rizzato 
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Fig.3f: Erichsen, Oliveira, and Rizzato 
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